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Abstract 

This paper is the second in a series of studies working towards constructing a 
realistic, evolving, non-potential coronal model for the solar magnetic carpet. In 
the present study, the interaction of two magnetic elements is considered. Our 
objectives are to study magnetic energy build up, storage and dissipation as a 
result of emergence, cancellation, and flyby of these magnetic elements. In the 
future these interactions will be the basic building blocks of more complicated 
simulations involving hundreds of elements. Each interaction is simulated in the 
presence of an overlying uniform magnetic field, which lies at various orientations 
with respect to the evolving magnetic elements. For these three small-scale 
interactions, the free energy stored in the field at the end of the simulation 
ranges from 0.2 — 2.1 x 10 26 ergs, while the total energy dissipated ranges from 
1.3 — 6.3 x 10 26 ergs. For all cases, a stronger overlying field results in higher 
energy storage and dissipation. For the cancellation and emergence simulations, 
motion perpendicular to the overlying field results in the highest values. For 
the flyby simulations, motion parallel to the overlying field gives the highest 
values. In all cases, the free energy built up is sufficient to explain small-scale 
phenomena such as X-ray bright points or nanoflares. In addition, if scaled for 
the correct number of magnetic elements for the volume considered, the energy 
continually dissipated provides a significant fraction of the quiet Sun coronal 
heating budget. 
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1. Introduction 



The magnetic carpet is the term used to describe the small-scale photospheric 
magnetic field of the quiet Sun ( Title and Schrijver, 1998[ ). Magnetic features 



in the magnetic carpet are generally categorised into three main classifications: 
ephemeral regions, network features and internetwork features. The nature of 



these small-scale features has been extensively studied by many authors (Harvey 



and Martin, 1973[ |Martin, 1988| |Martin, 1990[ |Wang et al., 1995| |Schrijver et 
al., 1997||Chae et al., 2001[ |Hagenaar, 2001||Parnell, 2002| |Hagenaar, Schrijver 



and Title, 2003 de Wijn et al 



2008; Parnell et al. 

_I l_ 



2009} [Thornton and Parnell, 



2011 ). The main properties of the magnetic features are summarised in Table 1 



Supergranular motions (Simon and Leighton, 1964 Hagenaar, Schrijver and 



Title, 1997 Paniveni et al., 2004[ ) provide the dominant flow pattern within the 
magnetic carpet. In addition to this, magnetic carpet features may evolve via 
the processes of emergence, cancellation, coalescence and fragmentation. These 
evolution processes, coupled with the photospheric flow pattern, mean that the 



quiet Sun photosphere is highly dynamic. Hagenaar, DeRosa and Schrijver ( 2008 ) 



determined the photospheric flux replacement time to be just 1-2 hours. Since 
magnetic fields from the magnetic carpet extend up through the chromosphere 
and into the lower corona, it is expected that the quiet Sun corona is also highly 
dynamic. 

Within the quiet Sun corona, energy may be released in a variety of ways. 
Withbroe and Noyes (1977) determined that the amount of heating required to 
maintain the quiet Sun corona is 3 x 10 5 ergs cm~ 2 s . There are a number of 
small-scale, transient phenomena associated with sporadic energy release, such as 
X-ray bright points (XBP) or nanoflares. XBPs are localised brightenings within 
the quiet sun corona, observed as point-like features or small loops (Golub et al., 



1974). They are typically associated with opposite polarity magnetic features on 



the photosphere, and in most cases with cancellation events (Webb et al., 1993). 



While XBP are not believed to be the primary so urce of quiet Sun heat ing, 

- 1 (iHabbal and Grace, 19911), 



as they radiate only around 5 x 10 4 ergs cm 2 s 



they do provide a contribution (Longcope and Kankelborg, 1999). Many authors 



have considered the interaction between pairs of small-scale magnetic elements in 



association with XBP (e.g. Priest, Parnell and Martin (1994), Longcope (1998), 



Longcope and Kankelborg (19991, von Rekowski, Parnell and Priest (2006)). 



Nanoflares are localised, impulsive bursts of energy, releasing energy on the order 
of 10 24 ergs (Parker, 1988). It is believed that nanoflares provide a contribution 



to coronal heating, by dissipating energy through magnetic reconnection. several 



authors have considered nanoflares in this context, including |Cargill| (|19 93), 

and 



Browning et al. 


Vekstein 


(2009 


)• 



(2008) 



Sakamoto, Tsuneta and 



In the past, several models for the solar magnetic carpet have been pro- 
duced. These fall roughly into two categories, although there is a natural overlap 
between them. The first are models that aim to reproduce the evolution of 



photospheric magnetic features as seen in magnetogram observations ( 


Schrijver 


et al., 1997 


van Ballegooijen et al., 1998 Simon, Title and Weiss, 2001 


Parnell, 


2001 Cranmer and van Ballegooijen, 2010 


Meyer et al., 2011 


). The second are 
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Table 1. Magnetic flux features within the solar magnetic carpet. 



Name 


Origin 


Flux (Mx) 


Size (km) 


Ephemeral 


Emerge within 


10 18 - 10 20 


3 000 - 5 000 


Regions (ERs) 


supergranules 






Internetwork 


Emerge within 


10 1G - 2 x 10 18 


2 000 (mean) 


(IN) Features 


supergranules 






Network 


Formed from residuals of 


10 18 - 10 19 


1 000 - 10 000 


Features 


ERs (90%) and IN features (10%) 







those that simulate the small-scale coronal field, usually using static photospheric 



boundary conditions ( Schrijver and Title, 2002) Close et al., 2003| Close et 
al., 2004; Regnie r7 Parnell and Haynes, 2008| |Cranmer and van Ballcgooijen, 



2010). These models tend to use the potential held approximation to study the 



coronal field. Previous modelling has shown that coronal heating inferred from 
the number of null points deduced from potential field models is insufficient to 



explain the hot corona (Schrijver and Title, 2002 Regnier, Parnell and Haynes, 



2008). As a result, more complex models that follow the history and distribution 



of electric currents and non-potentiality within the coronal field are required. 

With this in mind, our long term goal is to produce a non-potential model for 
the continuous evolution of the magnetic carpet, that models both the photo- 
spheric and coronal field. To construct this, we have formulated a two component 
model: a 2D photospheric model acts as the lower boundary for the full 3D coro- 
nal simulation. The photospheric component includes the effects of supergranular 
flows in addition to the four main flux evolution processes of emergence, cancella- 
tion, coalescence and fragmentation. In Meyer et al. (2011 ), hereafter Paper I, we 
show that our photospheric model reproduces many observed properties of the 
magnetic carpet. The purpose of the current paper is to present and test a new 
modelling technique that follows the continuous evolution of the 3D magnetic 
carpet coronal field through non-linear force-free states. This is in contrast to the 
models described above, which use static potential field extrapolations. A key 
difference between our non-potential model and the above potential field models 
is that our model retains a memory of flux connectivity and magnetic interac- 
tions over relevant timescales. We apply this technique to commonly occurring 
magnetic carpet interactions. These interactions are emergence, cancellation, and 
flyby. All three are simulated under the presence of an overlying magnetic field. 
Studying these interactions with the overlying field is useful, since small-scale 
phenomena on the solar photosphere can be influenced by stronger overlying 
fields. In a following paper, we will carry out full 3D magnetic carpet simulations 
of hundreds of magnetic elements, using the synthetic magnetograms described in 
Paper I as photospheric boundary data. Therefore the present paper is designed 
to consider three basic building blocks of the more complicated simulations. 
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Previous studies have also considered such interactions, these include full 



MHD simulations of magnetic flux emergence (e.g. Archontis et al. (20041, Mac- 
Taggart and Hood (2009)) and of a magnetic fiyby (e.g. Galsgaard, Parnell and 
Blaizot (20001, Parnell, Haynes and Galsgaard| (2010)). In the paper of Gals- 
gaard and Parnell ( 2005J), the authors study the heating associated with a fiyby 
using full MHD simulations, and they show that the amount of energy stored or 
dissipated within the corona depends on several factors. They conclude that it 
is not sufficient to know the evolution of the magnetic field at the photospheric 
boundary to be able to predict the evolution of the coronal field. Knowledge of 
the strength and direction of the overlying field is also important. 

Within the literature, there are no previous studies that compare all three 
interactions of emergence, cancellation and fiyby to one another under the same 
modelling approximations and using the same parameters. Therefore within the 
present study, we will consider the flux connectivity of the magnetic elements, 
the build up and storage of free magnetic energy, and the energy dissipated 
during these three basic magnetic carpet interactions. In particular, we wish to 
identify the factors involved in determining the amount of magnetic energy that 
is stored or dissipated. We use this study as a preliminary analysis to quantify 
results expected in more complex simulations involving more magnetic elements 
(see Paper I). 

The paper is outlined as follows. Section [2] discusses our treatment of the 
photospheric boundary condition and how it couples to the coronal magnetic 
field model. The set-up and analysis of the three basic interactions are described 
in Section [3] We present the discussion and conclusions in Section |4j 



2. Model 

A non-linear force-free magnetic field is a useful approximation to the coronal 
field, as it allows for the existence of electric currents and free magnetic energy. 
For a magnetic field to be force-free, it must satisfy: 



V x B = a(r)B, 



(1) 



and 



V • B = 0. 



(2) 



The physical conditions required for such an approximation to be valid are 



described by Regnier (2007). The parameter a(r) describes the twist of the 



magnetic field. It is a scalar function of position, but must be constant along a 
given magnetic field line. An important property of a non-linear force-free field 
is that it allows for regions of both high and low a, along with varying sign of 
a, so may therefore model a wide variety of coronal structures. 

There are several methods for constructing non-linear force-free fields from 
fixed photospheric boundary conditions; a summary is given by |Schrijver et 



al. (2006). In contrast to these methods, which produce single independent 



extrapolations, we choose to model a continuous evolution of the coronal field 
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through a series of quasi-static, non-linear force- free equilibria, that are driven by 
an evolving photospheric boundary condition. The magnetofrictional technique 
employed to carry this out is described next. 



2.1. Magnetofrictional Method 



Our model for the 3D coronal magnetic field is based upon that of van Ballegooi 



jen, Priest and Mackay (2000); the method is described fully in Mackay, Green 



and van Ballegooijen (20111. It has been used in the past to study the global 



coronal magnetic field ( Yeates, Mackay and van Ballegooijen, 2008), the evolving 



magnetic structure of solar filaments ( Mackay and van Ballegooijen, 2009 ) , and 



the decay of an active region (Mackay, Green and van Ballegooijen, 2011). In 



the model, we evolve the coronal field by the induction equation, 



dA 



= v x B + e, 



(3) 



where A is the vector potential and B — V x A. The plasma velocity v is 
specified by: 

lj xB 

V 

v B 2 ' 



(4) 



following the magnetofrictional method of Yang, Sturrock and Antiochos ( 1986 ), 
where v is the coefficient of friction and j = V x B. Through using this formal- 
isation, the coronal field evolves through a series of approximate, quasi-static, 
non-linear force-free equilibria. The second term on the right hand side of Equa- 
tion [3] is a non-ideal term that represents hyperdiffusion. The form is chosen to 
be: 



e = 



B 

£2 



(5) 



(Boozer, 1986 Bhattacharjee and Hameiri, 1986 Strauss, 1988), where 774, the 
hyperdiffusivity constant, is chosen to be 4.7 x 10 5 km 4 s _1 . The key effect of 
hyperdiffusion is that it conserves magnetic helicity whilst smoothing gradients 
in a. For any non-linear force-free field, hyperdiffusion acts to reduce the field 
towards a linear force-free state containing the same magnetic helicity. It should 
be noted that the timescales of the present simulations are far too short for a 
linear forcc-frcc field to be reached. 



2.2. Magnetic Energy Storage and Dissipation 



To consider the effect of the magnetic carpet on the corona, one aspect of the 
simulations that we are interested in is the build-up and release of energy. At 
any instant in time in our numerical box of volume V , the total magnetic energy 
is 

W = I v ^ dV - (6) 
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Following this, the rate of change of the total magnetic energy is 

1 f d ( B 2 * 



dW 



Substituting in Equation [3j 
dt{ 2 



47T 



B 



dt 



dB 

~dt 



dV. 



= B • V x (v x B + e) 



(7) 



where 



and 



dW 

~dT 



47T 



IdS - 



Q dV, 



(v x B + e) x B + r) A B 2 aSIa 



B 2 



Q = — Hv^+ 77 4 |Va|- 



(8) 



(9) 



(10) 



The first term on the right hand side of Equation [8] represents the energy 
injected due to surface motions, along with that injected or removed due to 
flux emergence or cancellation respectively. In addition, there is a contribution 
from hyperdiffusion. The second term is the rate at which energy is dissipated 
per unit time, due to the coronal evolution. This dissipated energy, which is 
released in the coronal volume, may be considered as energy that is available to 
be converted into heat or plasma motions. While it may be regarded as this, for 
simplicity in the present paper we only consider the size of Q, and do not follow 
the corresponding plasma processes. 



For the coronal dissipation term (Equation 10), the first term represents 
energy dissipation due to magnetofriction, which is released as the field relaxes 
to a force-free state. The second term represents energy dissipation due to hy- 



perdiffusion, which is described by van Ballegooijen and Cranmer (2008). In the 



present simulations, it is found that Va and the relaxation velocity are both 
largest near the magnetic elements and at locations where the magnetic field 
lines reconnect. Hence the dissipative term has its largest contribution at the 
sources, where there is strong magnetic field, B, and at locations of changing 
magnetic topology. 

Another quantity that we study is the free magnetic energy stored within the 
magnetic field. This is defined to be the difference between the energy of the 
non-linear force-free field and the corresponding potential field. The free energy 
at any instant is 



E f (t) = W Ql (t) - W p (t) 



dV, 



(11) 
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where B n ^ is the non-linear force-free field and Bp is the corresponding potential 
field. Within the simulations we will study first the energy that is continually 
released and may contribute to the heating of the corona (Equation [To]). Second, 
we will consider the energy stored in the magnetic field that may be related to 
more sporadic events, such as XBPs or nanoflares (Equation [TT]) . The relative 
importance of these two forms of energy will be discussed for each of the three 
interactions between magnetic elements. 

2.3. Photospheric Boundary Condition 

To simulate the processes of emergence, cancellation and flyby, each magnetic 
element at the photosphere is assumed to have a simple Gaussian form. Therefore 
the normal field component (B z ) of an element is given by: 

B z = B e- r2 / r o, (12) 

where B is the peak strength, r the Gaussian half-width and r the distance 
from the centre of the element. Each magnetic element within our simulations 
has a peak strength of Bq = 88 G, a Gaussian half-width of vq = 0.6 Mm and 
an absolute flux of 10 18 Mx. The total photospheric magnetic field is made up 
of the sum of a number of these elements. 

A unique feature of the simulations is that B z during the evolution of the 
photospheric field is specified analytically at discrete time intervals, Tj (= 200 
s = 3.3 min) apart. Movement of the sources between these time intervals is 
obtained by changing their central positions, (x, y), rather than advecting them 
numerically. Through doing so, we avoid undesirable numerical effects such as 
overshoot due to numerical differentiation, and pile-up at cancellation sites due 
to forcing. 

Using this description for magnetic elements, we may model a wide range of 
magnetic flux interactions, such as: 

(a) Cancellation: Two magnetic elements, one positive and one negative, move 
together until their centres coincide. If they are of the same strength, they 
completely cancel. A cancellation occurring over a time period of 100 min 
(30 Ti) can be seen in Figure [ij a). 

(b) Emergence: Two magnetic elements, one positive and one negative, have 
the same strength and initially the same central position. Moving the 
sources apart then simulates emergence as is seen in observations. An 
emergence occurring over a time period of 100 min (30 Tj) can be seen 
in Figure [ijb) . 

(c) Flyby: Two magnetic elements that move relative to one another, but never 

interact at the photospheric level. A flyby occurring over a time period of 
166.7 min (50 Tj) can be seen in Figure [l|c). 

At any instant in time, the flux through the photosphere is given by 

Flux= / B z dS, (13) 
Js 
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(a) 



(b) 



(i) t = 50 (ii) t = 66.7 (iii) t = 83.3 (iv) t = 96.7 (v) t = 100 



(c) 



(i) t = (ii) t = 3.3 (iii) t = 16.7 (iv) t = 33.3 (v) t = 50 



(i) t = (ii) t = 43.3 (iii) t = 83.3 (iv) t = 123.3 (v) t = 166.7 

Figure 1. Three magnetic flux interactions modelled at the photospheric level. The inter- 
actions occur between a positive and a negative magnetic element of equal strength: (a) 
cancellation, (b) emergence and (c) flyby. In each case, red and blue contours represent the 
positive and negative polarities respectively, at levels of ±[4,7,14,28,57] G. The area shown is 
of size 10 Mm X 10 Mm. For each image, the time is given in minutes. 




u- 1.0x10 



5.0x10 



40 60 
Time (min) 



100 



Figure 2. Total absolute flux through the photosphere as a function of time for a cancellation 
(black), emergence (blue) and flyby (red) event. 
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where S is the photospheric boundary surface. Figure [2] shows a plot of the 
total absolute flux through the photosphere as a function of time for each of the 
events in Figure [l] It can be seen that the curves are all completely smooth. The 
cancellation curve (black line) is level until the magnetic elements encounter one 
another (around t — 60 min); the curve then decreases smoothly to zero at which 
point the elements completely cancel. The emergence curve (blue line) shows the 
opposite behaviour, while the flyby curve (red line) remains completely level as 
no flux emergence or cancellation occur^] 

To ensure that V • B = within the simulations, the coronal field induction 
equation is specified in terms of A. Therefore to drive the evolution of the coronal 
field we require as photospheric boundary data the A x and A y corresponding to 
the analytically specified B z at each time step T{ . Without loss of generality, we 
can write A x and A y at z = in terms of a scalar potential $(x, y), where 

(44 V ) M =VX(^)=(- ,-— j. (14) 

Then from the z-component of B = V x A we have: 

<9 2 $ <9 2 $ . . 

B ^o) = -g^- W - (15) 

This may be solved using a fast Fourier transform to find and hence the A x 
and A y corresponding to B z at the level of the photosphere. 

To produce a continuous time sequence between each interval Tj, where B z 
and subsequently A x and A y are analytically specified, a linear interpolation 
of A x and A y is carried out between each Tj and T^+i, where 500 interpo- 
lation steps are taken between each analytical specification. This means that 
at each time Tj, the normal field at the photospheric boundary matches the 
exact analytically specified B z given by the sum of the Gaussian profiles of the 
discrete magnetic elements. In response to the photospheric evolution of A x and 
Ay at z=0, the vector potential within the coronal volume, and therefore the 
coronal field, evolves through a series of quasi-static equilibria as described by 
Equation [3j This treatment of the magnetic field at the photospheric boundary 
ensures that we still have freedom for A z , which sits half a grid point up from 
the photosphere (due to a staggered grid). Any non-potentiality in the coronal 
field near the photosphere arises from the z component of the vector potential. 
The initial condition for the coronal field in each simulation is a potential field 
extrapolated from A x and A y on z=0, at t = Tq. 

Since the evolution of our field is continuous, connectivity within the coronal 
magnetic field is maintained throughout the simulation. The resultant series 
of non-potential equilibria retain a memory of flux interactions from one step 
to the next, and subsequently the build up of non-potential effects. This is a 
significantly different approach compared to independent potential field extrap- 
olations, and provides a new insight into the energy budget of the quiet Sun 
corona. 



1 Notc that for the flyby, only 4=0-100 min is plotted, however the curve remains completely 
level throughout the entire simulation (4=0-166.7 min). 
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(a) 





(b) 



Figure 3. Cancellation: The red arrows represent the direction of the overlying field, the blue 
arrows represent the direction of motion of the magnetic elements. The bipole's axis is oriented 
(a) parallel to, (b) anti-parallel to, and (c) perpendicular to the overlying field. 



3. Basic Interactions of Magnetic Elements 

The three basic interactions studied - cancellation, emergence and flyby - are 
processes that commonly occur between magnetic elements of equal but opposite 
flux. Each simulation uses the photospheric boundary treatment described in 
Section |2.3| which is then applied to the 3D magnetofrictional model in order 
to drive the coronal field evolution. Section [3~T1 describes the features of the set- 



up that are common to all three cases; Sections 3.2 - 3.4 then give the results 



for individual cases. Section 3.5 compares the three basic interactions to one 
another. 



3.1. Set-up 

A numerical box of size 30 x 30 x 17.58 Mm is chosen, composed of 256 x 
256 x 150 grid cells. In each case, the interaction between magnetic elements is 
centred around the midpoint of the box. The box is periodic in the x-direction, 
and closed in the y-direction and at the top. 

For each interaction, the simulation is run under the presence of an overlying 
uniform magnetic field of strength 1 G, 5 G or 10 G, which points in the re- 
direction. Each interaction is also simulated with three different orientations of 
the bipole's axis with respect to the overlying field (parallel, anti-parallel or 
perpendicular). The effects of varying the strength of the overlying field and 
direction of motion of the bipole with respect to the overlying field are inves- 
tigated. For each simulation, the magnetic elements are advected at a constant 
velocity of 0.5 km s _1 . 

3.2. Cancellation 



The initial set-up for each of the cancellation simulations is illustrated in Fig- 
ure [3] Each magnetic element is initially positioned 3 Mm from box centre. It 
then takes 100 min to reach the midpoint, where the opposite polarity elements 
cancel. As seen in Figure [2j the total absolute flux through the photosphere 
(black line) remains constant until the magnetic elements come into contact, at 
which point the flux rapidly decreases to zero. 
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(c) (d) 

Figure 4. Cancellation of a bipole that is aligned parallel to a 5 G overlying field, (a) Non- 
linear force-free field and (b) potential field as seen in the x-y plane at z=0. (c) Non-linear 
force-free field and (d) potential field as seen in the x-z plane at y=15. A selection of magnetic 
field lines originating from the bipole at the photospheric level is plotted in each image. In 
images (c) and (d), some of the overlying field lines have also been plotted. The images are 
taken at t = 50 min. Red and green contours represent positive and negative magnetic field. 

3.2.1. Field Lines 

Typical examples of the magnetic connectivity during cancellation, for a 5 G 
overlying field and each orientation of the bipole can be seen in Figures [4] 
(parallel), [5] (anti-parallel) and [6] (perpendicular) , at t — 50 min. Images (a) 
and (c) show the non-linear force-free field. Although the field configurations for 
each case are only shown at one time, similar configurations occur up until the 
point at which the magnetic elements cancel. For parallel cancellation, no matter 
the strength of the overlying field, all flux from the positive polarity connects to 
the negative polarity. In contrast, for the anti-parallel and perpendicular cases, 
as the strength of the overlying field is increased, connectivity between the two 
magnetic elements decreases. 

In Figures |4j [5] and [6j images (b) and (d) show potential field extrapolations 
at t = 50 min, for the same photospheric field distribution as in images (a) 
and (c). For each case, magnetic field lines are plotted from the same starting 
points as in images (a) and (b). For the parallel cancellation (Figure [4]), the non- 
linear force-free field lines are similar to those of the potential field. However, for 
the anti-parallel and perpendicular cancellation, more twisting and bending of 
the magnetic field lines can be seen. This indicates that for the anti-parallel and 
perpendicular cases, the non-linear force-free field produces significantly different 
results. 
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(c) (d) 

Figure 5. Cancellation of a bipole that is aligned anti-parallel to a 5 G overlying field, (a) 
Non-linear force-free field and (b) potential field as seen in the x-y plane at z=0. (c) Non-linear 
force-free field and (d) potential field as seen in the x-z plane at y=15. A selection of magnetic 
field lines originating from the bipole at the photospheric level is plotted in each image. In 
images (c) and (d), some of the overlying field lines have also been plotted. The images are 
taken at t = 50 min. Red and green contours represent positive and negative magnetic field. 

3.2.2. Flux Connectivity and Energetics 

Figure [TJa) shows a plot of the total flux connecting the magnetic elements as a 
function of time, for parallel (black), anti-parallel (blue) and perpendicular (red) 
cancellation and a 5 G overlying field. Both the non-linear force-free field (solid 
lines) and corresponding potential field extrapolation (dashed lines) have been 
plotted for comparison. The total flux connecting the magnetic elements does 
not vary significantly between the non-linear force-free field and potential field 
cases, because the bipole is simply shrinking from an initially potential state, 
where many of the properties of the initial potential field are preserved. The 
most flux connects between the magnetic elements for the parallel case, and the 
least for the anti-parallel case. Similar results are found for the 1 G and 10 G 
overlying field simulations. However, as the overlying field strength increases, the 
amount of flux connecting the magnetic elements decreases in the anti-parallel 
and perpendicular cases. 

Figure ^f\b) shows the total magnetic energy as a function of time, for the 
perpendicular cancellation and a 5 G overlying field. The red solid line shows the 
non-linear force-free field energy, while the red dashed line shows the energy of 
the corresponding potential field. The non-linear force-free field energy initially 
increases as energy is injected due to surface motions, whereas the potential 
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(c) (d) 

Figure 6. Cancellation of a bipole that is aligned perpendicular to a 5 G overlying field, (a) 
Non-linear force-free field and (b) potential field as seen in the x-y plane at z=0. (c) Non-linear 
force-free field and (d) potential field as seen in the x-z plane at y=15. A selection of magnetic 
field lines originating from the bipole at the photospheric level is plotted in each image. In 
images (c) and (d), some of the overlying field lines have also been plotted. The images are 
taken at t = 50 min. Red and green contours represent positive and negative magnetic field. 



field energy is continually decreasing. However, both curves decrease as the 
magnetic elements begin to cancel, and there is an outflow of energy through 
the photospheric boundary. There is a slight increase in the non-linear force-free 
field energy towards the end of the simulation, as strongly curved overlying field 
lines form during the final stages of cancellation of the bipole. In contrast, the 
overlying field is completely straight in the final potential field. 

Figure[7|c) shows the variation of the free magnetic energy, -Ef(t), as a function 
of time, as computed using Equation |ll| Results are shown for the 5 G overlying 
field, and the line colours show the same orientations as in Figure [7]Ja). The 
free energy stored ranges from 0.37 — 1.57 x 10 26 ergs. The parallel cancellation 
results in the least free energy, while the perpendicular cancellation results in the 
most. On comparing Figures ^b) and (c), one can see that the free magnetic 
energy built up within each simulation is small compared to the total energy 
within the box (around 1% for the 5 G perpendicular case). However, in order to 
avoid boundary effects, the computational box is large compared to the magnetic 
elements and their area of evolution. Therefore the evolution of the magnetic 
elements only perturbs a small volume of the overlying field, and the majority 
of the overlying field remains in a near potential state. If we compare the free 
magnetic energy built up to the potential field energy contribution of the bipole, 
we see that the build up of free energy is significant. From Figure [Tib), we 
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Figure 7. Plots as a function of time for a cancelling bipole moving parallel (black), anti-par- 
allel (blue) and perpendicular (red) to a 5 G overlying field, for the non-linear force-free field 
(solid) and corresponding potential field (dashed), (a) Total flux connecting magnetic elements, 
(b) total magnetic energy, (c) free magnetic energy, Ef(t), (d) energy dissipated, f v QdV, and 
(e) cumulative energy dissipated, E^(t). (f) Free magnetic energy (triangles) and cumulative 
energy dissipated (stars) at the end of each simulation, for a 1 G (black), 5 G (blue) and 10 
G (red) overlying field. 



compute the potential energy contribution of the bipole to be 2.43 x 10 26 ergs. 
This is of similar order to the free energy stored in the 5 G perpendicular case. 

In addition to free magnetic energy stored within the system, energy is con- 
tinually being dissipated, as described by the heating term Q (Equation 10 1. 
This is illustrated by Figure ^d) , which shows Q integrated over the volume at 
a given instant in time. Again, results are shown for the parallel (black), anti- 
parallel (blue) and perpendicular (red) cancellation, with a 5 G overlying field. A 
greater quantity of energy is dissipated towards the end, when full cancellation 
occurs. More energy is dissipated in the anti-parallel and perpendicular cases 
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(c) (d) 

Figure 8. Cancellation of a bipole aligned anti-parallel to a 5 G overlying field: normalised 
contours of (a) Qf rc and (b) Q^d as seen in the x — z plane at y = 15 Mm, and normalised 
contours of (c) Qf rc and (d) Qhd as seen in the y — z plane at i = 15 Mm. In (a), (b) 
and (c), contour levels are [0.9,0.7,0.5,0.3,0.1]; in (d) contour levels are 200 times less 
[0.0045,0.0035,0.0025,0.0015,0.0005]. A selection of magnetic field lines originating from the 
bipole and the overlying field is plotted in each image. 



than in the parallel case. This is because smaller gradients in a are present in 
the parallel case as the bipole is simply shrinking and no reconnection occurs 
with the overlying field. 

The dissipated energy may be split into two terms representing energy dissi- 
pation due to magnetofriction, Qf rc , and hyperdiffusion, Qhd, where 

B 2 B 2 
Qfrc = -j— Hv| 2 and Q hd = — ?7 4 | Vct| 2 . 

Figure [8] shows normalised contours of Q {lc and Q hd in the x — z plane (Fig- 
ures [8j a) and (b)) and the y— z plane (Figures |8jc) and (d)) taken at t = 50 min 
for the anti-parallel cancellation with a 5 G overlying field. In each case, the plots 
are taken through the centre of the bipole. For (a) - (c), the contours are at the 
same levels, in (d) they are 200 times lower. Considering Figures |8ja) and (c), it 
can be seen that Q {rc occurs throughout a large part of the coronal volume, both 
along field lines connecting the magnetic elements, and along those connecting 
the magnetic elements to the overlying field. In general, the term is largest low 
down, where the field lines are perturbed by the motion of the magnetic elements. 
Also, in Figure (8^c), two 'wings' of Q frc can be seen suspended in the coronal 
volume. These originate due to strongly curved field lines that reconnect between 
the bipole and the overlying field. 

In contrast, from Figures|8jb) and (d), it can be seen that Q h d is more localised 
than Qf rc ; the reason is that Qhd only arises where gradients in a occur. This 
happens mainly at two locations: first, near the footpoints of the field lines 
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connecting the magnetic elements, and secondly around the separatrix surface 
which separates the overlying field lines from the closed connections between 
the magnetic elements. In all of the simulations discussed in this paper, we see 
a similar trend as to where energy dissipation occurs. In each case, Qf IC is seen 
to be space filling, whereas Q hd is more localised. 

By comparing Figures [7|c) and (d) , we see that J v QdV (expressed in ergs 
s _1 ) is around three orders of magnitude smaller than the free energy (expressed 
in ergs) stored by the end of the simulation. Although it is three orders of 
magnitude smaller, the values shown in Figure [7^d) are instantaneous values. 
By integrating Q over both the volume and time, we can see how much energy 
has been cumulatively dissipated over the whole simulation: 



E d (t) 



QdV 



v 



dt. (16) 



A plot of the total energy dissipated as a function of time is shown in Figure [7^e), 
for the 5 G case of each orientation (lines are coloured as in[7^a)). Cumulatively, a 
significant amount of energy has been dissipated. By comparing Figures[7|c) and 
(e) it can be seen that by the end of the simulation, more energy is cumulatively 
released (1.71 — 2.49 x 10 26 ergs) than is stored as free magnetic energy. 

Although Figures [7|a)-(e) only show results for the 5 G overlying field cases, 
Figure [7]jf) compares the values of free energy (triangles) and total energy 
dissipated (stars) at the end of the 1 G (black), 5 G (blue) and 10 G (red) 
simulations. The highest value of free energy is found for a 10 G overlying 
field, and perpendicular cancellation. In general, higher free energy is found 
for a stronger overlying field and when a greater volume of the overlying field is 
disturbed (i. e. perpendicular cancellation). We also see that a stronger overlying 
field results in more energy dissipation. For all overlying field strengths, parallel 
cancellation results in the least energy dissipation, however, a similar cumulative 
amount of energy dissipation occurs for both the anti-parallel and perpendicular 
cancellation. 

3.3. Emergence 

In the emergence simulations, initially, the net flux through the photosphere is 
zero. The magnetic elements, which coincide, then move apart until they each 
reach a separation distance of 3 Mm from the box midpoint after 100 min. As 
a result, a bipole appears in the photospheric distribution, simulating what can 
be classed in photospheric magnetograms as an emergence. The blue line in 
Figure [2] shows the total absolute flux through the photosphere as a function of 
time for the emerging bipole. Figure [9] illustrates the initial set-up for each of 
the simulations. 

3.3.1. Field Lines 



Figure 10 shows images of the 5 G emergence case, for each of the three ori- 



entations of the bipole's axis, at t = 50 min. At this time, the two magnetic 
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(b) (c) 

Figure 9. Emergence: The red arrows represent the direction of the overlying field, the blue 
arrows represent the direction of motion of the magnetic elements. The bipole's axis is oriented 
(a) parallel to, (b) anti-parallel to, and (c) perpendicular to the overlying field. 

elements have separated and no longer overlap. As in the cancellation simula- 
tions, when the emerging bipole's axis is parallel to the overlying field, all flux 
from the positive polarity connects to the negative polarity. In the anti-parallel 
and perpendicular cases, connections between the two magnetic elements exist in 
the early stages of emergence. However in both these cases, for strong overlying 
fields, connections between the magnetic elements can be completely severed by 
the end of the simulation. 

The photosphcric boundary conditions for the cancellation and emergence 
simulations are exactly the reverse of one another. This means that a potential 
field extrapolation at t = n min of each emergence case is identical to that at t = 
(100— n) min of the corresponding cancellation case. Thus, the same photospheric 
field distribution exists for both emergence and cancellation at t — 50 min. 
Therefore at this time, the field line plots of the non-linear force-free fields and 
potential fields in Figures |4{ [5] a nd [6] (cancellation) may be compared to the non- 
linear force-free fields in Figure |10| (emergence) . A comparison of these images 
shows that the emerging bipole's field is significantly different from that of the 
cancelling bipole or the corresponding potential field. In the parallel emergence 
case (|To]^c)), the low- lying field lines of the overlying field have been bent and 
pushed around either side of the bipole as the magnetic elements emerge. This 
shows that the special boundary treatment that we use allows the elements to 
emerge as a single flux system into the overlying coronal field. 

3.3.2. Flux Connectivity and Energetics 

Figure[TT|a) shows a plot of the total flux connecting the two magnetic elements 
as a function of time for parallel (black), anti-parallel (blue) and perpendicular 
(red) emergence with a 5 G overlying field. Both the non-linear force-free field 
(solid lines) and the corresponding potential field (dashed lines) are shown. For 
the perpendicular case, much more flux connects between the magnetic elements 
in the non-linear force-free simulations than in the potential field extrapolations. 
Due to the continuous nature of the magnetofrictional method, connections that 
are formed between the two magnetic elements as they first emerge are main- 
tained as they move apart. For the perpendicular case, the flux still connects 
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(f) 

Figure 10. Emergence of a bipole in the presence of a 5 G overlying field: (a), (b) and (e) 
as seen in the x-y plane at z=0, (c), (d) and (f) as seen in the x-z plane at y=15. For each 
case, the images are taken at t = 50 min. The bipole's axis is aligned: (a) and (c) parallel to, 
(b) and (d) anti-parallel to, and (e) and (f) perpendicular to the overlying field. A selection 
of magnetic field lines originating from the bipole at the photospheric level is plotted on each 
image. In images (c), (d) and (f), some of the overlying field lines have also been plotted. Red 
and green contours represent positive and negative magnetic field. 
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Figure 11. Plots for an emerging bipole moving parallel (black), anti-parallel (blue) and 
perpendicular (red) to a 5 G overlying field, as a function of time: (a) total flux connecting 
magnetic elements, (b) free magnetic energy, Ef(i), (c) energy dissipated, f v QdV , and (d) 
cumulative energy dissipated, E^it). 



between the magnetic elements by the end of the non-linear force-free field sim- 
ulation, but all connections have been completely severed in the potential field 
extrapolation. A comparison of the perpendicular simulation in Figure 11 ^a) 
to the corresponding cancellation plot (Figure [7^ a), red line) indicates that, 
for this orientation, the emergence shows a much larger departure from the 
potential field. Less departure is found for the anti-parallel case, since much more 
reconnection occurs as the bipole emerges into the oppositely aligned overlying 
field. 

Figure |ll[ b) shows the free magnetic energy as a function of time, for each 
orientation of the bipole with a 5 G overlying field (lines are coloured as in 
Figure [llja)). The plot shows that perpendicular emergence results in the most 
free energy, and the parallel emergence in the least, where energy values range 



from 0.32 — 1.12 x 10 26 ergs. Figure 11 'c) shows a plot of the rate of energy 



dissipation, Q, integrated over the volume as a function of time. In all cases, 
the energy dissipation rate increases rapidly as the bipole first emerges, and 
the photospheric flux is increasing. Subsequently, the rate of energy dissipation 
continues to increase as the two magnetic elements move apart, but the increase 
is less rapid. As in the cancellation simulations, if we compare Figures [iT|b) and 
(c), it can be seen that the instantaneous energy dissipation is three orders of 
magnitude less than the free energy by the end of the simulation. Figure [TT)(d ) 
shows the cumulative energy dissipated as a function of time (Equation 16), 
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Figure 12. Free magnetic energy (triangles) and cumulative energy dissipated (stars) at the 
end of each emergence simulation, for a 1 G (black), 5 G (blue) and 10 G (red) overlying field. 



where values range from 1.82 — 2.32 x 10 26 ergs. In contrast to Figure 11 'b) which 
shows that a perpendicular emergence results in the greatest build-up of free 
energy, Figures [TT^ c) and (d) show that more energy is dissipated per unit time 
and cumulatively in an anti-parallel emergence. This happens because that in 
the anti-parallel case, larger gradients in a are produced and more reconnection 
takes place as the bipole emerges into the oppositely aligned overlying field. 

Figure [12] compares the values of free energy (triangles) and total energy 
dissipated (stars) at the end of each simulation. The amount of free magnetic 
energy at the end of each emergence simulation follows the same trends as in 
the cancellation simulations: a stronger overlying field results in the build up 
of more free energy. The perpendicular emergence tends to result in the most 
free energy being built up, the parallel emergence in the least. We see that the 
most energy is dissipated in the anti-parallel case, and the least in the parallel 
case, as much more reconnection occurs in the anti-parallel case compared to 
the other two. As with the free energy, a stronger overlying field leads to more 
energy being dissipated. 

To test the results of our simulations, it is possible to estimate the maximum 
amount of free magnetic energy that can be built up and stored for the anti- 
parallel emergence simulations. If we assume that no reconnection occurs, a 



current sheet (see e.g. Parnell and Galsgaard (2004), Archontis, Tsinganos and 



Gontikakis (2010)) would separate the bipole's field from the oppositely directed 



overlying field. The free energy is then given by: 

F - B ^ D (17) 

where Bo is the strength of the overlying field, (j> is the absolute flux of each 
magnetic element and D is the photospheric separation of the magnetic elements. 



For the 5 G simulation, Equation 17 gives 4.8 x 10 ergs. This is the theoretical 



maximum value for the free energy built up. We note that within our simulations, 
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Figure 13. Flyby: The red arrows represent the direction of the overlying field, the blue 
arrows represent the direction of motion of the magnetic elements. The motion of the positive 
magnetic element is (a) parallel to, (b) anti-parallel to, and (c) perpendicular to the overlying 
field. 



the theoretical maximum cannot be obtained due to numerical diffusion and 
the fact that reconnection occurs as soon as the bipole emerges. We find that 
the occurrence of reconnection results in a smaller amount of free energy being 



stored. From Figure 12 the amount of free energy stored at the end of the 5 
G anti-parallel simulation is approximately 20% of the theoretical maximum. 
Although the free magnetic energy is one part of the energy calculated in the 
simulation, we also compute the energy dissipated due to relaxation processes 
and hyperdiffusion. For the 5 G case, a further 2.3 x 10 26 ergs of magnetic energy 
is cumulatively dissipated (Figure [Ti|d)). Summing the stored energy and the 
total energy dissipated, we obtain 3.2 x 10 26 ergs, which is close to the theoretical 
maximum of 4.8 x 10 26 ergs. Therefore, by taking into account not only the free 
energy, but also the energy dissipated when reconnection is allowed, we find a 
similar total amount energy to the theoretical maximum when reconnection is 
not allowed. 

3.4. Flyby 

The photospheric boundary distribution for the flyby simulations is slightly dif- 
ferent from the previous simulations. The flyby simulations are run for 166.7 min 
instead of 100 min, and each magnetic element is advected a distance of 5 Mm 
rather than 3 Mm. The total absolute flux through the photosphere is constant 
throughout the simulation (Figure [2j red line). The two magnetic elements are 
advected past one another under the presence of an overlying field, so that their 
final positions mirror their initial positions. An illustration of the initial set-up of 



each flyby simulation is shown in Figure 13 In two set-ups the magnetic elements 



are advected past one another in the x-direction. In the first case (Figure |l3[a)) 
the motion of the positive magnetic element is parallel to the overlying field, in 
the second case (Figure [l3^b)) it is anti-parallel to the overlying field. In the 
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(i) (ii) (iii) 

Figure 14. A series of x-y plane images for each flyby case with a 5 G overlying field. The red 
and green contours represent positive and negative magnetic field. The images in each case are 
taken at (i) t = min, (ii) t = 83.3 min and (iii) t = 166.7 min, and a selection of field lines 
is plotted originating from the magnetic elements. The positive magnetic element is advected 
(a) parallel to, (b) anti-parallel to and (c) perpendicular to the overlying field. 

third set of simulations (Figure [l3^c)) the magnetic elements are advected in 
the y-direction, so that their motion is perpendicular to the overlying field. 

3.4-1- Field Lines 

Figure [m] shows a series of x-y plane images at t = 0, 83.3 and 166.7 min for (a) 
parallel, (b) anti-parallel and (c) perpendicular flyby and a 5 G overlying field. 
Since the parallel flyby is the reverse of the anti-parallel flyby, the photospheric 
flux distribution at the start of each case is identical to that at the end of the 
other. This means that the initial potential fields in Figures [l4|^a)(i) and (b)(i) 
may be compared with the final non-linear force-free fields in Figures |l4|b) (iii) 
and (a) (iii) respectively. It is clear that the non- linear force- free field in each case 
is quite different from the corresponding potential field. In particular, significant 
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differences can be seen for the parallel flyby. For the non-linear force-free field 
( (a) (iii) ) , strong connections exist between the magnetic elements. However, 
there are no such connections in the potential field ((b)(i)). The field lines 
connecting the two elements are very twisted, and the bipole's magnetic field 
appears to occupy a much larger volume of the corona than in the corresponding 
potential field (Figure [l4^b)(i)). 

It is also of interest to compare Figures |T4^a) (ii) and (b)(ii), as their pho- 
tospheric distributions are also identical to one another. Even though both 
simulations have been running for the same amount of time (83.3 min), the shape 
of the bipole's field is very different. For the parallel simulation, the field lines 
that connect the magnetic elements are much more twisted. This illustrates that 
when the evolution of the coronal magnetic field is continuous, the properties of 
the field very much depend on its previous evolution and connectivity, not only 
on the photospheric boundary distribution. 

For the perpendicular flyby (Figure [l4{c)), the photospheric boundary dis- 
tributions in (c)(i) and (c) (iii) are symmetric to one another. At the midpoint, 
the bipole's axis becomes aligned fully parallel to the overlying field, and this 
results in two major occurrences of reconnection. The first occurs as the magnetic 
elements move towards one another, and the total flux connecting from one mag- 
netic element to the other rapidly increases. The second occurs after the magnetic 
elements have passed one another at the midline, causing connections between 
the magnetic elements to break, and the total connecting flux decreases. More 
flux remains connected between the two magnetic elements in Figure |14[ c) (iii) 
than in the corresponding potential field in Figurc[l4{c)(i). Again, this illustrates 
the effect of the memory of previous connectivity in our simulations. 

3.4-2. Flux Connectivity and Energetics 

Figure [l5|a) shows a plot of the total flux connecting the magnetic elements 
as a function of time, for a parallel (black), anti-parallel (blue) and perpen- 
dicular (red) flyby, with a 5 G overlying field. Both the non-linear force-free 
field (solid line) and corresponding potential field (dashed line) are shown. For 
the parallel case, at all times, the total flux connecting the magnetic elements 
is greater for the non-linear force-free field than the corresponding potential 
field. In particular, flux still connects between the magnetic elements at the end 
of the simulation, but does not in the corresponding potential field. For the 
anti-parallel case, initially, no flux connects between the magnetic elements. As 
the magnetic elements are advected past one another, the bipole's axis rotates 
to become aligned with the overlying field, and the total flux connecting the 
magnetic elements increases. Even though the connecting flux increases, the 
total flux connecting the elements is less than that of the potential field. 

For the potential field case of the perpendicular flyby (red dashed line), the 
plot of the total flux connecting the magnetic elements is symmetric about the 
line t = 83.3 min. In contrast, the connecting flux plot for the non-linear force- 
free field is non-symmetric. Before the polarities have aligned with the overlying 
field, the flux connecting the magnetic elements is lower than that of the potential 
field. However, after the elements pass one another at the midpoint, more flux 
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Figure 15. Plots for a flyby, where the positive polarity is advected parallel to (black), 
anti-parallel to (blue) and perpendicular to (red) a 5 G overlying field, as a function of time: 
(a) total flux connecting the magnetic elements, (b) free magnetic energy, E^(t), (c) energy 
dissipated, f v QdV, and (d) cumulative energy dissipated, E^(t). 



connects between them in the non-linear force-free field than in the potential 
field case. This once again indicates significant differences between the potential 
field extrapolations and the non-linear force-free field simulations which retain 
a memory of flux connectivity. 

Figure [l5| b) shows a plot of the free magnetic energy as a function of time for 
the 5 G case of each orientation (lines are coloured as in Figure [f5| a)), where the 
final values range from 1.50 — 2.05 x 10 26 ergs. The free energy stored by the end 
of the parallel flyby is greater than that of the anti-parallel flyby because, in the 
parallel case, more flux connects between the magnetic elements for a greater 
amount of time. However, the free energy is greatest for the perpendicular flyby, 
for two reasons. First, a greater volume of the magnetic field is disturbed by 
the magnetic elements, building up more free energy. Second, as the elements 
move past one another, flux from the positive magnetic element is forced to 
connect to the negative element. Numerous closed connections form and, due to 
the continuous nature of the magnetofrictional evolution, these flux connections 
are maintained as the magnetic elements move apart. Any free magnetic energy 
may be stored along them. 

Figure [l5|c) shows a plot of the rate of energy dissipation, Q integrated over 
the volume, as a function of time. In each case, the energy dissipated increases 
as the magnetic elements move towards one another, then levels off after they 
have passed one another at the midline. Initially, the curves for the anti-parallel 
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Figure 16. Free magnetic energy (triangles) and cumulative energy dissipated (stars) at the 
end of each flyby simulation, for a 1 G (black), 5 G (blue) and 10 G (red) overlying field. 



and perpendicular cases are very similar. However, towards the end of the simu- 
lations, the rate of energy dissipation becomes greater in the perpendicular case. 
This happens because in the perpendicular case the amount of flux connecting 



the magnetic elements is still changing. This can be seen in Figure 15 a) , where 
towards the end of the simulation, the slope of the total flux connecting the 
magnetic elements is much steeper in the perpendicular than in the anti-parallel 
simulation. Figure fl5[d) shows the cumulative energy dissipated as a function of 



time, as calculated by Equation 16 One can see from Figures |l5[c) and (d) that 
the energy dissipation is greatest in the case of the parallel flyby (5.67 x 10 26 
ergs in total) and least in the case of the anti-parallel flyby (4.70 x 10 26 ergs in 
total) . 

Figure [16] compares the values of free magnetic energy (triangles) and total 
energy dissipated (stars) at the end of each flyby simulation. The most free 
energy tends to be stored in the 5 G case of each orientation, with the greatest 
amount of free energy resulting from the 5 G perpendicular flyby. Between 2 — 3 
times more energy is cumulatively dissipated by the end of each simulation 
than is stored as free energy. The greatest amount of energy is dissipated for 
the parallel flyby with a 10 G overlying field. The exact amount of free energy 
stored and energy dissipated by the end of each simulation depends on a balance 
between the strength and orientation of the overlying field, the volume of the 
coronal field that is disturbed, the amount of rcconncction that occurs and, 
finally, the amount of connections that exist between the magnetic elements 
throughout their evolution. 

In the next section, we compare the cancellation, emergence and flyby simu- 
lations to one another. 

3.5. Comparison of interactions 

For each of the three bipole interactions, three different orientations of the inter- 
action with respect to an overlying magnetic field and three different strengths of 
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Figure 17. Free magnetic energy, Ef(t) (triangles), and cumulative energy dissipated, E^(t) 
(stars), at t = 100 min in each simulation, for a 1 G (black), 5 G (blue) and 10 G (red) overlying 
field, in the parallel (— >), anti-parallel (<— ) and perpendicular (_L) cases. The green dashed lines 
separate the plot into three columns representing the cancellation (left), emergence (middle) 
and flyby (right) simulations. 



overlying field have been considered. In this section, we compare the free energy 



and energy dissipated for all cases. Figure 17 shows a plot comparing the free 
magnetic energy (triangles) and total energy dissipated (stars) at t — 100 min for 
each simulation. The plot is split into three columns representing the cancellation 
(left), emergence (middle) and flyby (right) simulations. For each interaction, the 
parallel (—>), anti-parallel (<—) and perpendicular (_L) cases are shown with a 1 G 
(black), 5 G (blue) and 10 G (red) overlying field. Clearly, for the emergence and 
cancellation simulations, the smallest amount of free magnetic energy and energy 
dissipated arises in the parallel cases, while the anti-parallel and perpendicular 
cases lead to the largest energy stored and dissipated. In contrast, for the flyby, 
the parallel cases result in the most free energy and energy dissipated. However, 
for all cases the stronger the overlying field, the larger the energy values tend to 
be. 

The amount of free energy stored in the magnetic field at t = 100 min in each 
simulation varies between 0.2 — 1.9 x 10 26 ergs, where the greatest value arises 
for the perpendicular cancellation with a 10 G overlying field. As discussed in 



Section 3.2 these values of free energy are small compared to the total energy 
within the volume (between 8.7 x 10 26 — 6.3 x 10 28 ergs depending on overlying 
field strength). However, the box is large compared to the size of the bipole 
and the bipole's area of interaction. In reality on the Sun, many such magnetic 
elements would exist within such an area, with many of them continually inter- 
acting with one another. If we consider the power law distribution of |Parnell et| 



al. (2009), we can determine the number of small magnetic elements that would 



exist in the simulation region. Taking 
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where Nt = 3.6 x 1CP 17 cm -2 and </>o = 10 16 Mx, the number of magnetic 
elements expected in a 30 x 30 Mm region with flux in the range 10 17 — 10 18 Mx 
is 46, and the number with flux in the range 10 16 — 10 18 Mx is 374. Therefore in 
more complex simulations with a realistic number of magnetic elements, the free 
energy could be 1 — 2 orders of magnitude higher. Even though this is the case, 
for the simulations here, the amount of free energy built up in a single event is 
sufficient to explain sp oradic events such as nanofl ares, which release energy on 
the order of 10 24 ergs (Golub and Pasachoff, 1997). 

The total amount of energy dissipated by t = 100 min varies between 1.3 — 
3.2 x 10 26 ergs, where the simulation that results in the most energy dissipated 
is the parallel flyby. We may compute the average rate of energy dissipation in 
ergs cm~ 2 s _1 for each simulation. We consider the central region in each case 
(x = 10 — 20 Mm, y = 10 — 20 Mm) to focus in on the interaction location. We 
find that the average energy dissipation for each simulation occurs in the range 
1.50 — 4.95 x 10 4 ergs cm~ 2 s _1 . On comparing these values to the radiative losses 
of either t he quiet Sun (3 x 10 5 ergs cm " 2 s _1 ) or of a coronal hole (8 x 10 5 ergs 



Withbroe and Noyes ( 19771), we find that the rate of dissipation for a 



single event is too low. However, note that there would be many more magnetic 
elements on the Sun in a region of the size we have considered. Using the figures 
from above, 41 magnetic elements of flux 10 16 — 10 18 Mx would be expected, 
and this may increase the energy dissipation rate by an order of magnitude to 
observed levels. In addition, we may compare the radiative losses to those of an 
XB P. The radiative losses of a n XBP have been measured as 5 x 10 4 ergs cm~ 2 
s _1 ( jHabbal and Grace, 1991 ). This implies that the average energy dissipation 
in some of our simulations are of the correct size to explain the occurrence of 
such an event. 



4. Discussion and Conclusions 

The aim of this paper was to investigate the coronal consequences of three basic 
photospheric magnetic interactions: cancellation, emergence and flyby. Each in- 
teraction was simulated in the presence of an overlying uniform magnetic field, 
which was taken to be parallel, anti-parallel or perpendicular to the motion of the 
magnetic elements. The bipoles considered here were representative of small-scale 
photospheric magnetic features such as ephemeral regions or network features. 
Each bipole's physical extent was on the order of 3 000 — 4000 km and its total 
absolute flux was 2 x 10 18 Mx. 

In all cases, the 3D coronal magnetic field was initially in a potential state. A 
continuous evolution of the coronal field was then produced via a magnetofric- 
tional relaxation technique that evolved the field through a series of quasi-static, 
non-linear force-free equilibria in response to applied photospheric boundary 
motions. Our treatment of the photospheric boundary evolution was discussed 
in Section 12.31 The continuous nature of the coronal evolution means that con- 
nectivity within the coronal volume was maintained from one step to the next. 
In many cases, this continuity allowed connections to exist longer than those 
found in potential field extrapolations. This allowed free energy to be built up 
and stored along closed field lines. 
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For each of the simulations, two forms of energy were studied. First there 
was the free magnetic energy which was stored in the non-potential magnetic 
field. This energy may be regarded as that available for sporadic coronal events 
such as XBPs or nanoflares. Second, energy was continually being dissipated, as 



described by Equation 10 (Q). This may be considered as energy that is available 
to be converted to heat or plasma motions, although for simplicity, here we did 
not follow the corresponding plasma processes. With the formalisation used in 
the present paper, the dissipated energy arose from the relaxation process em- 
ployed, along with hyperdiffusion. We found that Q mainly arose low down near 
the magnetic elements where the magnetic field departed most from a potential 
state, and at sites of changing magnetic topology. 

The amount of free energy stored at t = 100 min in each simulation ranges 
from 0.2 — 1.9 x 10 26 ergs. The cumulative energy dissipated in each simulation 
after the same amount of time is greater than the free energy stored; for each 
simulation, anywhere from 1.3 — 3.2 x 10 26 ergs of energy has been dissipated 
after 100 min. The upper limits to both these values are higher if we consider 
the values at the end of the flyby simulations, which run for 166.7 min (2.1 x 
10 26 ergs for free energy, 6.3 x 10 26 for dissipated energy). For cancellation and 
emergence, the amounts of free and dissipated energy are smallest when the 
motion of the magnetic elements is parallel to the overlying field, and largest 
when it is perpendicular. In contrast, for flyby, the amounts of free and dissipated 
energy are greatest in the parallel case. In all cases, a stronger overlying field 
tends to lead to greater energy storage and dissipation. The simulation that 
results in the most free energy is the perpendicular cancellation with a 10 G 
overlying field, while the simulation that results in the most energy cumulatively 
dissipated is the parallel flyby with a 10 G overlying field. The exact amount of 
free energy stored and energy dissipated by the end of each simulation depends 
upon several factors: the strength and orientation of the overlying field, the 
volume of the overlying field that is disturbed, the amount of reconnection that 
occurs and, finally, the total flux connecting the magnetic elements. 

The free magnetic energy built up in the present simulations is small compared 
to the total magnetic energy within the volume. However, in each case, the free 
energy is a significant fraction of the bipole's energy contribution (8 — 86%), 
and is more than enough to account for small-scale, transient phenomena such 
as nanoflares or XBPs. The rate of energy dissipation in each case is between 
10 22 — 10 23 ergs s _1 . This could provide a contributio n towards the heating 
rate of an XBP of 3 x 10 23 — 10 24 ergs s _1 , determined by Habbal and Withbroe 



( 1981 ). We also find that for the inner 10 x 10 Mm of each simulation, the energy 
dissipation rate is between 1.50 — 4.95 x 10 4 ergs cm~ 2 s _1 . This is equivalent 
to 5 — 17% of the energy required to heat the quiet Sun corona (3 x 10 5 ergs 



cm ' s " 1 , 



Withbroe and Noyes (19771). Although it is at most 17% of the 



coronal energy budget, on the Sun, many tens or hundreds of such small-scale 
magnetic elements would be found in a region of the size modelled here. The 
continual interaction of these magnetic elements with one another would result 
in a significantly larger build up of free magnetic energy and greater energy 
dissipation. With the expected number of magnetic elements in a region of this 
size, the free energy and energy dissipation rate may easily be 1 — 2 orders of 
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magnitude larger than those found for the simulations in this paper. This would 
bring it in line with coronal requirements. It is therefore of key importance to 
consider more complicated simulations of multiple magnetic elements. 

From this study, the next step is to simulate the coronal evolution of the 
synthetic magnetograms constructed in Paper I. This would allow us to study 
the energetics of many events at the same time, as occur on the Sun. This will 
be the aim of the third paper in our series. In these more complex simulations, 
we will study many aspects of the coronal evolution such as locations of coronal 
null points, electric currents, free magnetic energy and energy dissipation, and 
relate these to the dynamic processes occurring in the photospheric evolution. It 
will also be of interest to conduct a similar study by applying the same technique 
to real magnetogram data, such as those from SDO/HMI. Regions of interest in 
the resultant coronal magnetic field could then be compared with a wide variety 
of SDO coronal wavelengths. 
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